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The application of higher-order Laguerre Gauss (LG) modes in large-scale gravitational wave 
detectors has recently been proposed. In comparison to the fundamental mode, some higher-order 
Laguerre Gauss modes can significantly reduce the contribution of coating Brownian noise. Using 
frequency domain simulations we give a detailed analysis of the longitudinal and angular control 
signals derived with a LG33 mode in comparison to the fundamental TEMqo mode. The performance 
regarding interferometric sensing and control of the LG33 mode is found to be similar, if not even 
better in all aspects of interest. In addition, we evaluate the sensitivity gain of the implementation 
of LG33 modes into the Advanced Virgo instrument. Our analysis shows that the application of the 
LG33 mode results in a broadband improvement of the Advanced Virgo sensitivity, increasing the 
potential detection rate of binary neutron star inspirals by a factor 2.1. 
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I. INTRODUCTION 

The search for gravitational waves (GW) has led to 
a new class of extremely sensitive laser interferometers. 
The first generation of large-scale laser-interferometric 
gravitational wave detectors [TJ |5J |3J 2] is now in opera- 
tion with the aim of accomplishing the first direct detec- 
tion of gravitational waves. The detector performance is 
limited by several fundamental and technical noises. In a 
constant effort the noise contributions are minimised to 
improve the detectors signal-to-noise ratio. One of the 
limiting noise sources of the currently planned second 
generation gravitational wave detectors will be thermal 
noise [5 of the mirror test masses. There exist several 
components to thermal noise of which the Brownian ther- 
mal noise is largest in current interferometer topologies 
utilising arm cavities. Cooling of the mirror test masses 
as currently studied in CLIO [B] reduces the thermal noise 
provided an appropriate material is chosen for the optics. 
A different way to lower the thermal noise is to change the 
mode shape of the laser beam inside the interferometer. 
All current detectors use the fundamental TEMqo mode, 
but several other mode shapes such as Mesa beams ff], 
conical modes ^Sj and higher-order Laguerre Gauss (LG) 
modes [S] have been proposed for reducing thermal noise. 
The basic idea is to reduce thermal noise by generating a 
more uniform light intensity distribution on the mirrors 
without introducing higher clipping losses [9] . 

The proposed candidates for different beam shapes 
can be divided in two groups. The first group, which 
comprises flat-top and conical beams, would require the 
use of non spherical mirror shapes. As a result these 
modes are not compatible with current GW detectors 
and their spherical mirrors. Currently it is not clear to 
what precision these non-spherical mirrors can be man- 
ufactured and little experience in using such mirrors has 
been gained so far. The second group consists of higher- 
order Laguerre Gauss modes, which are fully compati- 



ble with spherical mirrors as the currently used TEMqo 
mode[lT]. So far LG modes have been mainly employed 
in the field of cold atom and quantum optics for example 
as optical tweezers [TD] or waveguides [11 . 

Currently several techniques for the generation of 
higher-order LG-modes exist e.g. using holograms [12l 
112], gratings [T3] and mode transformers [T3|TB]. With 
these techniques a conversion efficiency of 60% fTi has 
been demonstrated. Recently the creation of higher- 
order LG-modes with a very high mode purity [17 has 
been reported. Our paper assumes that using these tech- 
niques, higher-order LG-modes can be created with high 
power output and high mode purity required in the field 
of GW detection. We analyse the compatibility of such 
higher-order LG-modes with the core interferometer in 
future GW detectors, using Advanced Virgo in partic- 
ular as an example for a second-generation GW detec- 
tor. In Section |TT] we introduce the definition of the LG- 
modes, how they can be described in a Hermite Gauss 
mode basis system and how the coating brownian ther- 
mal noise is calculated depending on which LG-mode is 
used. Moreover we introduce some practical considera- 
tions concerning clipping loss, beam sizes and radius of 
curvature (RoC) of the mirrors which are essential for the 
later analysis. In Section [lTT| we perform a phase coupling 
analysis of a single arm cavity and a Michelson interfer- 
ometer using higher-order LG modes in comparison to 
the currently used fundamental mode. We determine any 
differences in their phase coupling between the different 
longitudinal and alignment degrees of freedom. In Sec- 
tion IIVI a numerical model based on a set of Advanced 
Virgo design parameters is used to analyse the prospects 
of higher-order LG modes in comparison to the currently 
proposed use of the fundamental mode. The detector 
sensitivities of the different interferometer configurations 
are computed to derive the envisioned detector inspiral 
ranges. In total we analyse and compare three different 
cases with each other. 
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II. INITIAL CONSIDERATIONS 



A. Hermite Gauss and Laguerre Gauss modes 

The Hermite Gauss (HQ) modes and Laguerre Gauss 
modes both present complete basis sets such that each 
LG mode can be presented by a sum of HG modes and 
vice versa. The so-called helical Laguerre Gauss modes 
can be written as [T51 [12] : 

Up,i (r,<^,z)= ^^^^^exp(i(2p+|;| + l)^'(z)) 



I / 2r^ 



expf-ifcsfe+i^^) (1) 



with r, and z as the cylindrical coordinates around 
the optical axis, the beam radius, \E'(z) the Gouy 

phase, q{z) the Gaussian beam parameter and L^^ix) the 
associated Laguerre polynomials. The indices must obey 
the following relations: < |/| < p where -p is the radial 
mode index and / the azimuthal mode index. 

The decomposition of these modes into Hermite Gauss 
modes can be performed as follows [20] : 

(2p+l) 

Mp,i(a;,2/,z) = ^ i''6(/ + p,p,fc)ug,';;_fcfe(a;,y,z) (2) 
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with real coefficients 
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where P^'^{x) denotes the Jacobi Polynomials. It is in- 
teresting to note in Equation l2| that a given LGpi mode 
is constructed of 2p -I- Z -I- 1 HG„„i modes of the order 
n + m = 2p + I. For example the LG33 mode is con- 
structed out of ten Hermite Gauss modes of the order 



B. Coating brownian thermal noise of 
Laguerre-Gauss modes 

According to [H] the power spectral density of dis- 
placement equivalent thermal noise is generally given by 



(4) 



with (f) being the loss angle and U the strain energy of the 
static pressure profile on the mirror surface normalised 
to 1 N. The interested reader is referred to [22j where 
detailed calculations of various thermal noises such as 
substrate brownian, coating brownian and thermoelastic 
thermal noise are presented. The currently most limiting 
thermal noise in GW detectors is the coating brownian 
thermal noise. In the case of a semi-infinite mirror the 



coating brownian thermal noise induced by a LGpj mode 
can be calculated using the strain energy 



p, Locating 
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Here 5c is the thickness of the coating, a is the Poisson 
ratio, Y is the Young modulus, w is the beam width at 
the mirror and gp^i is a scaling factor depending on the 
used LGpi mode. In the case of the fundamental LGqo 
mode this scaling factor is (70,0 = 1? whereas for a LG33 
mode 53,3 = 0.14 has to be used. Hence the power spec- 
tral density of displacement equivalent coating brownian 
thermal noise is more than a factor of seven smaller for 
a LG33 mode in comparison to the fundamental LGqo 
mode. For finite mirror sizes [32] finds that the deviation 
to the semi-infinite case are very small if the clipping loss 
of the beam on the mirror surface are small. 



C. Clipping loss and beam scaling factors 

For our later analysis it is important to know the clip- 
ping loss Zciip that affects the propagating Gaussian mode 
at the mirror due to its finite size. It is given by 

lciip{w,p,z) = 1-1 d(j) dr-r 
Jo Jo 
X u{w, r, (j), z) u* {w, r, <p, z) , (6) 

were w is the beam radius at the mirror, p is the radius 
of the mirror coating and u{w, r, 0, z) is the transversal 
field distribution of the mode of interest. Please note 
that the parameter beam radius w is a, measure of the 
beam size of the fundamental Gaussian mode (LGqo or 
HGoo). Higher-order LG or HG modes of the same beam 
radius actually are more spatially extended, in the sense 
that a significant amount light power can be detected 
at distances off the optical axis larger than the beam 
radius. Using the general definitions of the transversal 
field distribution for Hermite Gauss and Laguerre Gauss 
modes u we can compute the clipping losses of any of 
these modes. 

In Figure [1] the clipping losses for the fundamental and 
two higher-order LG modes are plotted over the mirror- 
radius to beam-radius ratio. One can see clearly that in 
comparison to the fundamental LGoo mode, the higher- 
order LG modes have a much more widely spread inten- 
sity distribution for a given beam radius. Hence they 
require either larger mirrors, or reduced beam radii for a 
fixed mirror size. Table |T] comprises the respective scal- 
ing factors, normalised to a optimised mirror size for a 
LGoo^niode with a clipping loss of 1 ppm. 
One consequence of this more widely spread intensity dis- 
tribution of the higher-order LG modes is of major impor- 
tance for the later analysis: in order to fit a higher-order 
mode optimally on the same mirror as the fundamental 
mode, the beam radius of the higher-order mode must 
be different from that of the fundamental mode. This 
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LGoo 


LG33 


LG55 


mirror to beam radius ratio 


2.63 


4.31 


5.05 


relative mirror radius 


1 


1.64 


1.92 


relative beam radius 


1 


0.61 


0.52 



Table I: Comparison of the mirror to beam radius ratio for 
1 ppm clipping loss and the corresponding scaling factors for 
the beam radius and mirror radius to keep the clipping loss 
constant normalised to the LGoo-mode with a clipping loss 
of 1 ppm. 




Mirror to beam size ratio 



Figure 1: The plot shows the relative power loss in reflection 
of a finite sized mirror due to clipping loss for three different 
incident transversal modes over the mirror to beam radius 
ratio. The different curves shown are: blue curve = TEMoo 
mode; black curve = LG33 mode; red curve = LG55 mode. 
The mirror to beam radius ratio for a fixed clipping loss of 
1 ppm is given in Table |l] 

corresponds to a different wave front curvature and con- 
sequently to a different spherical curvature of the cavity 
mirrors |18j . Thus changing an existing optical experi- 
ment such as an interferometer from a configuration us- 
ing e.g. the TEMoo mode to a configuration using the 
LG33 mode, the radii of curvature of the mirrors must 
be changed, if one wants to keep the clipping losses at a 
constant level. In most cases this would necessitate in a 
complete exchange of the mirrors. 



III. PHASE COUPLING COMPARISON OF 
THE LG33 MODE WITH THE FUNDAMENTAL 
HGoo MODE 

The sensitivity of future gravitational wave detectors 
will be limited partly by thermal noise. The use of higher- 
order LG modes represents a very interesting option for 
reducing this limit. For a successful implementation, 
however, higher-order LG modes must comply with the 
stringent phase noise requirements in these detectors. In 
the following we compare the phase noise performance of 



the currently used HGoo niode with that of a LG33 mode 
which serves as a representative of the family of higher- 
order LG modes. The analysis was performed with the 
numerical interferometer simulation Finesse [2T , which 
uses the Hermite Gauss modal expansion for describing 
the spatial properties of light fields transverse to the op- 
tical axis. In order to simulate higher-order LG modes 
we used the decomposition of higher-order LG into HG 
modes presented in Section |ll] The results of the anal- 
ysis described in this section are in principle applicable 
to many other Laguerre Gauss modes of interest. For in- 
stance for the LG55 mode, the corresponding mirror and 
beam radius scaling values are given in table |l] 



A. Configurations of interest 

Our phase noise coupling analysis uses the currently 
planned 3 km long Advanced Virgo interferometer as 
testbed. We compare the use of a LG33 mode with two 
different configurations using a fundamental mode: 

1. ) The LGq™^'' configuration: This configuration uses 

the optical parameters presented in j24j and rep- 
resents our reference configuration. The configu- 
ration uses arm cavity mirrors with RoCs Rc — 
±1910 m which corresponds to a beam size of the 
fundamental LGoo niode of w; = 35.2 mm at the 
mirrors surface and a corresponding waist size of 
Wo = 16.3 mm. According to table [l] this config- 
uration would have a clipping loss of 1 ppm for a 
mirror radius of /Osmaii = 92.5 mm. 

2. ) The LG33 configuration: This configuration uses 

the higher-order Laguerre Gauss mode LG33. It 
shares its beam parameters with the reference 
LGqo"*'' configuration to simplify the comparison. 
The mirror radius has to be adapted for this con- 
figuration because of the more wider intensity pro- 
file of the LG33 mode. A mirror radius of piargc — 
151.8 mm is required to maintain a clipping loss of 
Ippm (see Table |l|. 

3. ) The LGoo'^° configuration: The third configuration 

uses the fundamental LGoo mode in combination 
with the larger mirrors with radius piargc used in the 
LG33 configuration. As a result the beam size on 
the mirrors can be increased to w — 57.7 mm while 
still maintaining a clipping loss of 1 ppm. Hence 
the waist size decreases to wq = 8.9 mm and all 
other beam parameters change accordingly. 

The beam parameters of each configuration are displayed 
in Table [n] and their transversal intensity distribution on 
the mirror surface is shown in Figure [2j It is worth not- 
ing, that a comparison between the LG33 and the LGoo"^^" 
configuration is much more reasonable because these two 
configurations use the same mirror sizes. 
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Figure 2: Transversal intensity distribution at the mirror surface for the three configurations under investigation. 





T /^small 


LG33 


T large 


Rc [m] 


1910 


1910 


1536.7 


w [mm] 


35.2 


35.2 


57.7 


Wo [mm] 


16.3 


16.3 


8.9 



Table II: Beam and mirror parameters of the three different 
configurations used in the analysis. There exist always two 
RoC settings for achieving a given spot size on the mirrors. 
We have chosen the cavity geometry which reduces the radi- 
ation pressure induced alignment instabilities |25] . 



B. Tilt to longitudinal phase coupling of a single 
cavity 

The first part of our analysis is performed for a single 
cavity. At first we want to investigate the longitudinal 
error signal of each configuration in order to find out 
how they compare against each other. To generate this 
longitudinal error signal we use the Pound-Drever-Hall 
technique based on a modulation/demodulation scheme 
[25] . We found that all resulting error signals are iden- 
tical, which confirms that such error signals only depend 
on the average phase of the beam, and is independent of 
its modal distribution. 

Next we analyse the coupling of misalignment of a cav- 
ity into longitudinal phase noise for the three configura- 
tions of interest. The first results are comprised in Fig- 
ure [3] which shows the intra cavity power over the tilt 
angle (3 of the end mirror (EMX) and the longitudinal 
tuning (j) of the cavity for each configuration. The tuning 
value (j) of the cavity is the result of a modulo division 
of the cavity length by the wavelength. In the following 
this tuning is given in degrees with 360° referring to one 
wavelength: 



360° • {L mod A) 



(7) 



In the analysis the tilt angle (3 of the cavity's end mirror 
— shown on the x-axis — was varied between /zrad to 



1 ^rad. On the y-axis the tuning (j) was chosen such that 
the resonance of the cavity is clearly visible. One can see 
that a tilt of the end mirror changes the tuning of the res- 
onance for every configuration. Because the tuning refers 
to the length of the cavity, there is indeed a coupling from 
tilt into the longitudinal phase degree of freedom of the 
cavity. The coupling strength is different for each con- 
figuration. Nevertheless, configuration LGq™^'' and LG33 
behave in a very similar manner. Both show a shift of 
the resonance condition of approximately = 0.4° for 
a tilt of 1 /irad. In contrast, configuration LGqq'^®'' shows 
an increased coupling strength by more than one order 
of magnitude. For a tilt of 1 /^rad the cavity tuning of for 
the LGqq'^^° configuration changes by A0 « 5°. 

The tilt of the end mirror (EMX) also changes the ge- 
ometry of the eigenmode of the cavity. Compared to the 
input beam the eigenmode of the cavity is tilted as well. 
As a result the mode matching efficiency into the cav- 
ity is decreased, leading to a reduced intra cavity power. 
This behaviour is hardly visible in Figure [3] but clearly 
shown in Figure [J] The LGq™''" configuration is the most 
robust in terms of decreased intra cavity power, followed 
by the LG33 and the LGqq'^^'' configurations. 

In conclusion of this section the LGp™ ^" configuration 
performs best in all aspects of the analysis. Nevertheless 
as stated initially a comparison between the other two 
configurations is much fairer because they share the same 
mirror size. Taking this into account, the favourable 
mode is the LG33, because its tilt induced coupling into 
the longitudinal phase and into the intra cavity power is 
much less than for a LGqo mode on the same mirrors. 



C. Alignment analysis of a single arm cavity 

The proper sensing and control of the alignment de- 
grees of freedom of a gravitational wave detector is crit- 
ical for its successful operation. Hence, a comparison 
analysis of the alignment error signals of the individual 
arm cavities for the three different transversal mode con- 
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Figure 3: Intra cavity power over tilt angle (5 of the end mirror (EMX) and longitudinal tuning (/> of a single cavity shown 
for all three configurations of interest. The first two configurations show almost the same coupling from tilt into longitudinal 
tuning which is more than an order of magnitude lower than the coupling of the third configuration. 
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Figure 4: Comparison between all three cases and their cou- 
pling between tilt of mirror EMX to the tuning of mirror 
IMX. 



diode current is then demodulated at frequency f2. The 
demodulation phase is chosen to maximise the error sig- 
nal slope for the corresponding mirror. An example of 
the corresponding alignment error signals for the mis- 
alignment of both arm cavity mirrors using the LGg™^'' 
configuration can be seen in Figure |6] A good way of 
comparing different alignment error signals is looking at 
the resulting control matrix [28, .29^ . In our example the 
control matrix contains the values of the alignment er- 
ror signal slopes a in the working point generated by the 
two quadrant photodiodes QPDl and QPD2 for a mis- 
alignment of the arm cavity mirrors IMX and EMX, see 
Figure |5]@. The subscript of a indicates the readout 
diode and the superscript refers to the mirror of inter- 
est. Hence, the control matrix is given in general by the 
following expression 



configuration 



^IMX ^EMX 
MX 



'QPDl 
^IMX 



(8) 



figurations defined in Section III A is needed and pre- 
sented in this section. Full alignment control systems of 
advanced GW detectors are very complex and depend on 
the details of the detector design. The concept however, 
is firmly based on the control of resonant cavities and 
we can use a single Fabry Perot cavity to test whether 
LG modes are compatible with current alignment control 
systems. The analysis uses an optical layout as shown in 
Figure [5]®. An electro-optic modulator (EOM) imprint- 
ing a phase modulation with frequency f2 in combination 
with two quadrant photodiodes is used to generate align- 
ment error signals for the two arm cavity mirrors using 
the Ward technique described in [27]. Here each quad- 
rant photodiode is responsible for obtaining an alignment 
error signal of one arm cavity mirror. In the following we 
only consider mirror rotations around vertical axis. The 
results however remain applicable for the tilt degree of 
freedom of each mirror. The optimisation of the error 
signals does not use the theoretical optimal parameters 
but is done by tuning the parameters; this reflects realis- 
tic experimental procedures. The two quadrant photodi- 
odes are placed such that their Gouy phase difference is 
90° leaving the total Gouy phase arbitrary. Each photo- 



Applying the optimisation procedure for the alignment 
signals described earlier for all three configurations of 
interest, results in the following three control matrices. 



LG3 



LG" 



5.6152 0.0477 
2.1607 3.5878 



5.6152 



1 



0.009 



0.385 0.639 



7.444 0.022 
2.741 4.771 



7.444 



1 



0.003 



0.368 0.641 



17.774 15.330 
11.472 2.725 



17.774 



1 0.862 
0.645 0.153 



(9) 



(10) 



(11) 



For an ideal control matrix all matrix elements on the 
diagonal would be one and the off-diagonal elements are 
zero. Comparing the resulting control matrices, one can 
see that the LGqq'^° configuration performs worse than 
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Figure 5: Different optical layouts used in the alignment analysis described in SectionjlirC] Left @: The generation of alignment 
error signals for a single arm cavity. Right (S): Michelson interferometer with differentially misaligned arm cavities to study 
the power coupling into the output port of the interferometer. 



the other two: both mirrors couple much more strongly 
into QPDl than into QPD2 making the alignment er- 
ror signals far from ideal. The LGg™''" and the LG33 
configuration show a much better and almost even per- 
formance. This is represented by the fact that in both 
of these configurations the misalignment of mirror IMX 
couples three times more strongly into QPDl than into 
QPD2. Any misalignment of mirror EMX couples a fac- 
tor of 75 stronger into QPD2 compared to the signal 
sensed with QPDl for the LGq™^" configuration. This 
factor increases further to 216 in the case of the LG33 
configuration. 

In conclusion we can say that the LGqq'^° is outper- 
formed by the other two. The reason for this is not to 
be found in the mode shape but in the cavity geome- 
try. The RoCs of the mirrors of the LG^^^" and the 
LG33 configurations are the same, which results in al- 
most the same control matrix. In contrast the RoCs of 
the mirrors of the LGqq'^° configuration are much smaller. 
For the LGgo'^° configuration one obtains a cavity g- 
parameterflS] oi g — 0.91 which is very close to the insta- 
bility border of unity. The other two configurations have 
a g-factor of 0.33 which corresponds to a much more sta- 
ble and robust geometry. 



D. Coupling of differential arm cavity 
misalignment into the output port power 

An important measure for the performance of the dif- 
ferent optical modes in a full interferometer configuration 
is how much the differential misalignment of the arm 
cavities and the corresponding position change of each 
cavity's eigenmode couples with the power in the inter- 
ferometers output port due to the reduce mode overlap 
on the beam splitter. A sketch of this behaviour can 
be seen in Figure |5](b). This coupling mechanism can 



generate a signal which is indistinguishable from a GW 
signal. Hence it is important to analyse how higher-order 
Laguerre Gauss modes compete with the currently used 
fundamental mode. The following analysis compares this 
coupling mechanism for the three configurations of in- 
terest with respect to a reference value. This reference 
value is the limit Advanced LIGO fSIT specifies for the 
differential arm length deviation which is supposed to be 
smaller than 10~^^m [21]. This deviation together with 
the envisaged dark fringe offset of 10""'^^ m |32] allows us 
to calculate the expected increase in power in the output 
port of the interferometer to be 7.124 • 10^^ W. This ref- 
erence value can now be compared to the output power 
enhancement of the three configurations of interest for 
a misalignment of up to 200 prad (see Figure [t]) . The 
LGq™'''' configuration shows the smallest coupling from 
differential tilt of the arm cavities into the output port 
power followed by the LG33 and finally the LGqq'^®'' con- 
figuration. The values where each configuration crosses 
the reference limit of 7.124 • 10~® W are summarised in 
Table ITTTl 



LG33 

T /-I large 



misalignment angle /3 



148 prad 
46 prad 
26 prad 



coupling scaling factor 



1 

3.2 
5.7 



Table III: Results of the coupling analysis between the differ- 
ential arm cavity misalignment and the interferometers out- 
put port power. For each of the three mode configurations 
of interest the misalignment angle is given at which the in- 
terferometer output power has risen from it nominal value 
of 6.412 • 10"^ W to the reference limit of 7.124 ■ 10"'' W. 
The given coupling scaling factor shows the relative coupling 
strength of each configuration referenced to the LGqo'^" con- 
figuration. 



In conclusion the LG^ 



configuration performs at 



7 




-1 -0.5 0.5 1 

IMX misalignment angle p [urad] 




-1 -0.5 0.5 

EMX misalignment angle |3 [urad] 



Figure 6: Comparison of the obtained alignment error signals 
in the LGq™"" case sensed with the two quadrant photodi- 
odes while mirror IMX (top) and mirror EMX (bottom) is 
misaligned. 



least a factor of 3.2 better than the other two configura- 
tions. Nevertheless if we compare the two configurations 
which use the same mirror size, we find that the LG33 
performes much better than the LGq"^° configuration. 

Overall we can conclude from the performed phase 
coupling comparison analysis that higher-order Laguerre 
Gauss modes are suitable for the implementation in fu- 
ture GW detectors. The direct comparison between the 
two configurations which use the same mirror sizes - the 
LGq'J^'®'' and the LG33 configuration - clearly shows that 
in all of the presented phase coupling analyses the con- 
figuration using higher-order Laguerre Gauss modes per- 
forms better than the configuration based on the tradi- 
tionally used fundamental mode. This serves to underline 
the great potential and prospects of using higher-order 
Laguerre Gauss modes in future GW detectors. 
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Figure 7: Comparison of the dark port power when the arm 
cavities are differentially misaligned (see Fig.[5](S)) in reference 
to the dark fringe power that results from a differential arm 
length deviation of 10^^^ m while the interferometers dark 
fringe offset is 10"^'^ m. 

IV. PROSPECTS FOR ADVANCED VIRGO 

In this section we focus on the currently planned sec- 
ond generation GW detector Advanced Virgo [33, and 
its sensitivity. The design efforts for Advanced Virgo are 
rapidly progressing, yielding frequently improved detec- 
tor configurations. Hence the configuration presented in 
the following is unlikely to be the final Advanced Virgo 
configuration, but rather represents a snapshot of a devel- 
opment process. The numerical computations of the de- 
tector sensitivity have been performed with GWINC 

m- 

In the following we compare the expected sensitivity 
for a configuration using the fundamental TEMqo mode 
against a configuration using a LG33 mode in three differ- 
ent scenarios in which we assumed a fixed mirror radius 
of 'Tmirror = 0.17m[331. The three scenarios are: 

1 . The beams of both mode configurations experience 
an identical clipping loss at the ITM/ETM mirrors. 
Hence, the beam size used in each configuration is 
different which results in a different RoC of the arm 
cavity mirrors, (see Sec. [h]). This means that for 
a change of the mode shape used in the interfer- 
ometer from e.g. initially the TEMqo mode to the 
LG33 mode, all interferometer mirrors have to be 
exchanged. 

2. Both mode configurations use arm cavity mir- 
rors with an identical RoC. This enables a simple 
switching between the two different mode shapes. 

3. The planned thermal compensation system (TCS) 
is used to change the RoC of the arm cavity mir- 
rors by a certain amount. Hence one can start with 
one mode configuration which uses arm cavity mir- 
rors with the optimal RoC. Later the TCS system 
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enables us to change the RoC to optimise the clip- 
ping loss of the other mode configuration. Hence, 
this scenario is divided into two parts. The first 
one starts with an optimised parameter set for the 
TEMqo mode and the TCS system is used to im- 
plement the possible LG33 mode configuration. In 
the second part of this scenario we will start with 
an optimised parameter set for the LG33 mode and 
than change to a TEMqo mode configuration using 
the TCS system. 

The figures of merit of the comparison are the resulting 
effective detection range for a binary neutron star inspiral 
Tns/ns ^iid effective detection range for a binary black 
hole inspiral rgfj/BH- 

A. Scenario (i): Identical clipping loss 

The two mode configurations which are compared in 
this scenario will be the reference configurations through- 
out the whole analysis because of their fixed clipping loss 
of 1 ppm. According to Section |TT] the beam sizes of each 
configuration are different as well as the RoC of the in- 
terferometer mirrors. This scenario and its two configu- 
rations allow us to see the potential of using higher-order 
Laguerre Gauss modes in future GW detectors for a fixed 
mirror size. 

Table |IV| comprises the input parameters used in the 
analysis as well as the resulting inspiral ranges. As one 
can see the analysis was performed for two different sig- 
nal recycling (SR) detunings [35] - 750Hz and 300Hz - 
to emphasise that the improvements vary with the de- 
tuning. For both detunings the improvement using the 
LG33 mode is significant. The total mirror thermal noise 
is decreased by a factor of up to 1.68. According to our 
simulation this results in a relative improvement of the 
inspiral ranges by 20%-21% for a SR detuning of 750 Hz 
and 23%-28% for a SR detuning of 300 Hz. Hence the 
potential event rate[ll] of the Advanced Virgo detector 
can be increased by a factor of up to 2.1 by using the 
LG33 mode instead of the fundamental TEMqo niode. 



B. Scenario (ii): Identical RoCs 

Our second scenario uses arm cavity mirrors with a 
fixed RoC leading to the same beam size at the mir- 
rors for both mode configurations. The advantage of 
this configuration is that it allows using either TEMqo 
or the LG33 mode without exchanging the interferome- 
ter mirrors. To achieve this with a reasonable sensitivity 
for both configurations, a tradeoff concerning the clip- 
ping loss has to be made, see Section |TT] for more details. 
As a result the beam size of the TEMqq mode will be 
smaller compared to the reference configuration. On the 
one hand this results in much smaller clipping loss for this 
configurations, but on the other hand the clipping loss of 



the LG33 mode configuration will go down to an accept- 
able value as well. The input parameters and the results 
of the simulation analysis for this scenario are shown in 
Table |V] The comparison of these configurations shows 
that the inspiral ranges in the case of the LG33 mode 
configuration are greater by at least 35% and up to 53%. 
Although this looks promising, these two configurations 
have to be compared with the reference configurations 
for the same SR detuning analysed in scenario (i). We 
then find that the LG33 configuration of scenario (ii) per- 
forms 0.1% better concerning the BH/BH inspiral range 
but the NS/NS inspiral range is 11% worse because of its 
higher clipping loss of Idip = 30 ppm due to the slightly 
larger beam size. This results in a decreased intra-cavity 
power which lowers the sensitivity of the detector. The 
TEMqo configuration of scenario(ii) is also much less sen- 
sitive compared to the reference configurations because 
of the small beam size and the resulting higher thermal 
noise. The inspiral ranges of the TEMqq configuration 
decrease by approximately 22% in comparison to the ref- 
erence configuration. The weak performance of the both 
configurations renders this scenario not very favourable 
to be implemented into Advanced Virgo. 



C. Scenario (iii): Use TCS to adapt RoCs 

Our third scenario combines the advantages of the two 
earlier scenarios - high performance and compatibility 
between the different mode configurations. We propose 
to use the TCS, which is an essential part of future GW 
detectors, to introduce a constant offset onto the RoC 
of the arm cavity mirrors. This technique has already 
been demonstrated in the GEO 600 detector to match 
the RoC of the two interferometer end mirrors [35]. The 
basic idea of this approach is to start with one mode con- 
figuration with a clipping loss of /clip = 1 ppm at each 
arm cavity mirror. The TCS will then be used later 
to change the RoC of the arm cavity mirrors to opti- 
mise the useable beam width for the corresponding other 
transversal mode. The currently planned TCS |37j for 
second generation GW detectors are designed to intro- 
duce a RoC change of ARc = -120 m. The TCS uses a 
ring heater placed near the mirror substrate. The ther- 
mal radiation produced, is partly absorbed by the mirror 
substrate thus deforming its original shape. By placing 
the ring heater either in front or behind the mirror one 
can actually change the sign of the RoC change allowing 
an adjustment of Ai?c ± 120m. The two different signs 
of the possible RoC offset allow two different approaches 
whose input parameters and results are comprised in Ta- 
ble IvD 

The first approach uses the reference TEMqq mode 
configuration with /clip = 1 ppni as initial interferome- 
ter configuration. This configuration with an arm cavity 
RoC of Rc = 1522.8 m has already been analysed in sce- 
nario (i) . A constant change of the arm cavity RoC by 
ARc = -1-120 m introduced by the TCS will minimise the 
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Scenario 


SR det. 


w in 


^clip 


RoC 


Tns/ns 


Tbh/bh 


(i) 


[Hz] 


[cm] 


[ppm] 


[m] 


[Mpc] 


[Mpc] 


TEMoo 


750 


6.47 


1 


1522.8 


139.83 


1135.2 


LG33 


750 


3.94 


1 


1708.4 


168.34 


1373.9 


TEMoo 


300 


6.47 


1 


1522.8 


148.85 


1076.2 


LG33 


300 


3.94 


1 


1708.4 


191.26 


1322.4 



Table IV: Input parameters and results of the GWINC simulation analysis of the scenario (i) which uses identical clipping loss 
at the arm cavity mirrors for each mode configuration. 



Scenario 


SR det. 


w in 


^clip 


Rc 


Tns/ns 


Tbh/bh 


(ii) 


[Hz] 


[cm] 


[ppm] 


[m] 


[Mpc] 


[Mpc] 


TEMoo 


750 


4.22 


8.O0-9 


1647.2 


110.86 


900.47 


LG33 


750 


4.22 


30 


1647.2 


149.25 


1375.3 



Table V: Input parameters and results of the GWINC simulation analysis of the scenario (ii) which uses mirrors with an 
identical RoC for each mode configuration. 



beam size of a LG33 mode. This minimised mode still ex- 
periences a clipping loss of IcUp = 40.9 ppm. Please note 
that in this case the clipping losses of the LG33 mode 
will go down further with a larger RoC offset. The major 
advantage of this approach is that we can use the refer- 
ence TEMqo mode configuration initially and than later 
exchange it for a better performing LG33 mode configu- 
ration. This is reflected in the resulting inspiral ranges of 
these two configurations (see Table VI I . We find that the 
NS/NS and the BH/BH inspiral ranges of the LG33 mode 
configuration are increased by 2% and 18%, compared to 
the reference TEMqo mode configuration. For a higher 
RoC change ARc we can expect a much better perfor- 
mance due to the lower clipping loss. To reach the ref- 
erence LG33 mode configuration from scenario (i) which 
has a clipping loss of l^np — 1 ppm, requires for example 
an induced RoC change by the TCS of ARc = 186 m. 

The second approach uses the TCS to start with the 
reference LG33 mode configuration with ^ciip — 1 ppm 
and an arm cavity mirror RoC of Rc = 1708.4 m (see sce- 
nario (i) ) . The TCS can be used to introduce a constant 
offset in the arm cavity mirror RoC of ARc = —120 m, 
which maximises the possible beam size of a correspond- 
ing TEMqo mode configuration to w — 4.71 cm corre- 
sponding to the very small clipping loss Idip =4.8e-6 ppm. 
If the TCS is able to change the RoC by a larger amount, 
the usable beam size of the TEMqo mode configuration 
could be increased further. This approach leads to an 
increase of the inspiral ranges of the LG33 mode con- 
figuration by between 42% and 43% compared to the 
TEMqo mode configuration. Despite these large improve- 
ments with the LG33 mode configuration it is interesting 
to compare the TEMqq mode configuration used here to 
the reference TEMqq mode configuration used in scenario 
(i) to see how much one would lose by using this second 
approach of scenario (iii). It turns out that the inspiral 
ranges go down by approximately 15%. 



In conclusion scenario (iii) is the most promising one 
and the idea to use the TCS to introduce a constant 
change to the RoC of the arm cavity mirrors has great 
potential. It would allow a change from a TEMqq mode 
configuration to a LG33 mode configuration without ex- 
changing the main optics of the interferometer. To decide 
which of the two approaches described is the better one, 
one has to judge if one could either afford to have clipping 
losses for the LG33 mode configuration of Idip = 40.9 ppm 
in the first approach, or to have a decreased performance 
of 15% by the TEMqq mode configuration in the second 
approach compared to the reference TEMqq mode con- 
figuration. 



D. Sensitivity improvement from the use of LG33 
on the example of Advanced Virgo 

As an illustrating example we evaluate how much the 
sensitivity of an detector like Advanced Virgo could be 
improved by the application of LG33 modes. In par- 
ticular we compare configurations featuring TEMqq and 
LG33 modes with identical clipping losses, but RoCs op- 
timised for the individual modes. Our analysis made use 
of the same parameters as the most recent and compre- 
hensive noise model for Advanced Virgo [31] , (35] . We as- 
sumed the following reduction factors for the LG33 mode 
in our GWINC simulation: The coating Brownian and 
the substrate Brownian noise are reduced by factors 1.7 
and 1.9^3], respectively, while the thermo-elastic noise 
increases by a factor 1.7 [40|. 

Figure [s] displays the resulting sensitivity curves (blue 
traces) of Advanced Virgo featuring TEMqq (dashed line) 
and LG33 (solid line) modes. As coating Brownian noise 
is directly limiting the Advanced Virgo sensitivity in the 
frequency range between 40 and 200 Hz, an impressive 
sensitivity improvement can be achieved by application 
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Scenario 


SR det. 


w in 


^clip 


Rc 


Tns/ns 


Tbh/bh 


(iii) 


[Hz] 


[cm] 


[ppm] 


[m] 


[Mpc] 


[Mpc] 


TEMoo 


750 


6.47 


1 


1522.8 


139.83 


1135.2 


LG33 


750 


4.25 


40.9 


1642.2 


142.93 


1345.6 


LG33 


750 


3.94 


1 


1708.4 


168.34 


1373.9 


TEMoo 


750 


4.71 


4.8e-6 


1588.2 


118.40 


960.98 



Table VI: Input parameters and results of the GWINC simulation analysis of the scenario (iii) which uses the same clipping 
loss at the arm cavity mirrors for each mode configuration. 
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LG33: Coating Brownian noise 
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control matrices for a single Fabry-Perot cavity as well as 
the coupling of differential arm cavity misalignment into 
dark port power for a full Michelson interferometer with 
arm cavities. We were able to show that the LG33 mode 
performs similar if not even better than the commonly 
used TEMoo for all considered aspects of interferometric 
sensing. This strongly indicates that all currently avail- 
able experience and technology for interferometric sens- 
ing and control, which is based on the TEMqo mode, can 
be transferred to the use of the LG33 mode. Changing 
over from the fundamental mode to the LG33 will not 
require any fundamental changes of the interferometric 
control strategy or the control hardware, but only small 
adjustments of the involved parameters, such as servo 
gains. 



Figure 8: Sensitivity improvement from the implementation 
of the LG33 mode into the Advanced Virgo detector. This 
analysis is based on the detector configuration described in 
the Advanced Virgo Preliminary design [38j . For simplicity 
only the contributions of coating Brownian noise (b, c) and 
quantum noise (a) are shown, while all other noise contribu- 
tions are omitted in the plot, but taken into account for the 
overall sensitivity (d, e). The sensitivity improvement when 
going from TEMqo (d) to LG33 (e) mode corresponds to an 
improvement of the binary neutron star inspiral range from 
145 to 185 Mpc and increase the detectors NS/NS event rate 
by a factor of 2.1. 



of the LG33 mode. This broadband sensitivity improve- 
ment is concentrated on the range from 30 to 400 Hz with 
a maximal gain around 75 Hz. The binary neutron star 
inspiral range increases from 145 to 185 Mpc, which cor- 
responds to a rise of the expected NS/NS event rate by 
a factor of 2.1. 



V. SUMMARY AND OUTLOOK 

We carried out a comprehensive analysis of the 
prospects of higher-order LG modes for future gravita- 
tional wave detectors. Using numerical interferometer 
simulations, we compared the behaviour of the LG33 
mode with the fundamental mode (TEMoo). Our anal- 
ysis included tilt to longitudinal phase coupling, gener- 
ation of angular control signals and the corresponding 



In addition, we performed a quantitative evaluation of 
the expected sensitivity improvement from application 
of the LG33 mode, using the planned Advanced Virgo 
detector as an example. Three different options of how 
to change over from the TEMoo to the LG33 mode have 
been developed. The first scenario considers replacing 
the main mirrors, by ones with radii of curvature op- 
timised for the LG33 mode, while the second scenario 
assumes that the same mirrors are used for both modes, 
resulting in different clipping losses. In the third sce- 
nario the Advanced Virgo thermal compensation system 
is used to adjust the mirror curvatures for the two optical 
modes of interest independently. The maximum sensitiv- 
ity improvement is found to be achievable when replacing 
the mirrors (scenario (i)). Using the latest design pa- 
rameters of Advanced Virgo we were able to show that 
the application of the LG33 mode can give a broadband 
improvement of the Advanced Virgo sensitivity for all 
frequencies in the range from 30 to 400 Hz. The cor- 
responding binary neutron star inspiral range increases 
from 145 to 185 Mpc, enhancing the potential detection 
rate of binary neutron star inspirals by a factor 2.1. 

The next steps towards the actual implementation of 
higher-order LG modes in future GW detectors have to 
include the demonstration of efficient generation of high 
power LG beams, followed by setting up of tabletop ex- 
periments for experimental verification of the simulations 
presented in the first half of this article. 
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